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We study statically homogeneous Bose-Einstein condensates with spatially inhomogeneous inter- 
actions and outline experimental realizations of compensating linear and nonlinear potentials that 
can yield such constant-density solutions. We illustrate how the presence of a step in the nonlinearity 
coefficient can be revealed dynamically by a defect-dragging experiment due to the inhomogeneity of 
the sound speed. We then use effective-potential theory to perform a detailed analytical investiga- 
tion of the existence and stability of solitary waves in this setting, and we corroborate these results 
computationally using a Bogoliubov-de Gennes linear stability analysis. We find that dark solitary 
waves are unstable for all step widths, whereas bright solitary waves can become stable through a 
symmetry-breaking bifurcation as one varies the step width. We illustrate the various scenarios that 
permit this bifurcation by examining the phase plane of the corresponding dynamics. 



Introduction. For more than two decades, Bose- 
Einstein condensates (BECs) have been a fruitful experi- 
mental, computational, and theoretical testbed for inves- 
tigating nonlinear phenomena. In the mean-field limit, 
a BEC is governed by the Gross-Pitaevskii (GP) equa- 
tion [T], which is a nonlinear Schrodinger (NLS) equa- 
tion with an external potential. The NLS is important 
in many fields [2], and many ideas from disciplines such 
as nonlinear optics have proven important for investiga- 
tions of BECs. The ability to control various parameters 
in the GP equation makes it possible to create a wide 
range of nonlinear excitations, and phenomena such as 
bright [3l|4], dark [SHZ], and gap 8] solitary waves (and 
their multi-component [9] and higher-dimensional |10[lll] 
generalizations) have been studied in great detail using 
a variety of external potentials PJJ [UJ . 

The GP equation's cubic nonlinearity arises from a 
BECs interatomic interactions, which are characterized 
by the s-wave scattering length. The sign and magni- 
tude of such interactions can be controlled using Fesh- 
bach resonances [121414"] , and this has led to a wide array 
of theoretical and experimental possibilities [3l 01 [T5l [16] . 
In a recent example, Feshbach resonances were used to 
induce spatial inhomogeneities in the scattering length in 
Yb BECs [17] . Such collisional inhomogeneities, which 
amount to placing the BEC in a nonlinear potential 
in addition to the usual linear potential, can give rise 
to effects that are absent in spatially uniform conden- 
sates [T8H20] . This includes adiabatic compression of 
matter waves [31], enhancement of the transmission of 



matter waves through barriers |22j . dynamical trapping 
of solitary waves |22j , derealization transitions of matter 
waves [23], and many other phenomena. 

In the present work, we study the situation that arises 
when spatial inhomogeneities in nonlinear and linear po- 
tentials are tailored in such a way that they compensate 
each other to yield a constant-density solution of the GP 
equation. We demonstrate how this scenario can be en- 
gineered in experiments and investigate it for a step-like 
configuration. This situation is particularly interesting 
because the inhomogeneity is not mirrored in the BECs 
density profile, which makes the step undetectable using 
static density measurements. We show, however, that 
the step can be detected dynamically using a dragging- 
impurity experiment [24) . and we observe the emission 
of dark solitary waves when the dragging speed is above 
a critical velocity (which is different inside and outside 
the step). Additionally, we use effective-potential the- 
ory to examine the existence and dynamical robustness 
of dark and bright quasi- ID solitary waves for various 
step-potential parameters. 

Model and Setup. Starting with the 3D time- 
dependent GP equation, we consider a cigar-shaped con- 
densate and average over the transverse directions to ob- 
tain a quasi-lD GP equation [JJ [TJJ [IT]. In perform- 
ing the averaging, we assume that the BEC is strongly 
confined in the two transverse directions with a trap- 
ping frequency of uj± [38]. The solution of the quasi- 
lD GP equation is a time-dependent macroscopic wave- 
function ^(z,i). We assume the standing- wave ansatz 



2 



Step 



No step 




40 60 80 100 
t 



20 40 60 80 100 
t 



FIG. 1: [Color online] Numerical computations of defect drag- 
ging in the quasi-lD GP equation. (Left) Emission of a 
dark solitary wave as a defect is dragged through the step. 
(Right) The same computational experiment but without a 
step (so there is no solitary- wave emission). The defect speed 
is v — 0.6, and the other parameter values are 7 = — 1 and 
AV = 0.5. 



^(z,t) = 4>(z)e lflt to obtain the time-independent GP 
equation 



<t> Z z- n0 + V cxt (z)(f> + g{z)\cj>\ 2 4 = O, (1) 



where (j) is measured in units of (2|ao|) -1 ^ 2 and 
g(z) is a spatially varying nonlinearity associated with 
the (rescaled) scattering length a(z) through g(z) = 
a(z)/\ao\- We measure length in units of a± = 
\Jhj (mwi) and time in units of wj 1 , where m is the 
mass of the atomic species forming the condensate. The 
constant is the value of the scattering length in 
the collisionally homogenous system. Equation ([!]) has 
two conserved quantities: the number of atoms N — 
(oj./[2|ao|]) \^\ 2 dz and the Hamiltonian [TT] . 

For a square-step linear potential, one can use the 
Thomas-Fermi approximation (cf> zz = 0) for the ground 
state [TT]. Equating the densities inside and outside of 
the step then gives the constraint 
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where Vq and go are the constant background linear and 
nonlinear potentials, and AV and Ag are the differences 
between the step and background values of V(z) and 
g{z). The parameter 7 thus measures (and balances) the 
relative strengths of the steps in the linear and nonlinear 
potentials. To preserve smoothness, we implement the 
steps using hyperbolic tangent functions: 

AV 

V(z) = V + AV(z) = V Q + — [tanh(z+) - tanh(z_)] , 

g( z )=go + Ag(z) = .9o + ^ [tanh(z + ) - tanh(z_)] , (3) 

where z± = [z ± zq)/s, the step width is 2zq, and s 
controls the sharpness of the step edges. From ([2]), it 
follows that AV — jAg. For the remainder of this Letter, 
we take V = and \gi\ — — 1. This yields 7 = — 1 
and corresponds to nonlinear and linear steps of equal 
and opposite depths/heights. 



Experimentally, a spatially varying magnetic field 
B(z) results in a proportionally varying linear poten- 
tial V(z) — mFgF^BB(z) for magnetic spin states (with 
magnetic quantum number rrip, Lande factor gp, and 
Bohr magneton /i^) at sufficiently low magnetic fields 
within the regime of validity of the linear Zeeman ef- 
fect. For specific atomic species and spin states, there 
is an additional resonance dependence (a Feshbach reso- 
nance) of the coupling constant g on the magnetic field: 

g(B) = gbg — b-b ) ' wnere 9bg is the background 
coupling constant, Bq is the resonance field, and A is the 
resonance width. The condition of compensating linear 
and nonlinear potentials is fulfilled within the Thomas- 
Fermi approximation when n = — ^ . In theory, this 
implies that for any given density n, there is a field B c 
near the resonance Bq matching this requirement. Conse- 
quently, the density must remain constant for any static 
profile B(z) as long as B(z) s» B c and g(B) is an ap- 
proximately linear function of B. In practice, however, 
large nonlinearities lead to fast three-body recombination 
losses from traps and hence have to be avoided An 
atomic species with appropriate properties is cesium, for 
which the above conditions are fulfilled at typical densi- 
ties of 10 13 — 10 14 cm -3 for fields near the narrow Fesh- 
bach resonances at 19.8 G and 53.5 G [29\. Optical dipole 
traps near the surface of atom chips provide an environ- 
ment in which magnetic fields can be accurately tuned to 
and varied about the critical magnetic fields B c at these 
values. 

Step Detection. We detect the presence of steps in the 
potentials by performing computational experiments in- 
volving the dragging of a defect through the BEC [24] , l2"6 ] . 
For the above linear and nonlinear steps, the condensate 
density is constant within and outside of the step. How- 
ever, the speed of sound c is different in the two regions: 
c = y/g(z)n(z), where n(z) = \(f)(z)\ 2 is the BEC density 
[2 7) . We simulate a moving defect computationally using 
a potential of the form V(z,t) — Ae^ z ~ r ^ l w , where 
A and w are amplitude- and width-related constants, and 
r(t) = r(0) + vt represents the center of the defect mov- 
ing with speed v. The dynamics of defects moving in 
a BEC are sensitive to the speed of the defect relative 
to the speed of sound: speeds in excess of the speed of 
sound (i.e., supercritical defects) lead to the formation of 
dark solitary waves travelling behind the defect, whereas 
speeds below the speed of sound (i.e., subcritical defects) 
do not [26] . 

When the speed is subcritical, there is a density depres- 
sion with essentially the same functional form as the lin- 
ear potential. This changes shape slightly in the presence 
of the step; it deepens and widens for a step with Ag < 0, 
and it becomes shallower and narrower when Ag > [39j . 
When the defect speed is supercritical within the step re- 
gion but subcritical outside, we expect the nucleation of 
dark solitary waves in the step region. Because the de- 
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feet's speed is smaller than the background sound speed, 
the emission of solitary waves downstream of the defect 
becomes a clear indication of the presence of a step. We 
demonstrate this second scenario in Fig. [T] The third 
possible scenario involves a defect that is supercritical in 
both regions. 

Existence and Stability Criteria for Solitary Waves. 
Our detection scheme for the step in the potential sug- 
gests the ability to spontaneously produce solitary mat- 
ter waves. We now examine the existence and stability 
of such solutions as a function of the step parameters 
(especially its width). 

We apply the Bogoliubov-de Gennes (BdG) ansatz 



i ^ {z)+Y J {n j {z)e- i ^+v*{z)e i ^)] (4) 
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to the time-dependent quasi-lD GP equation. Equation 
Q defines the linear eigenfrequencies ujj for small per- 
turbations characterized by eigenvectors Uj(z) and Vj(z). 
Linearizing the time-dependent GP equation about the 
reference state <j>o{z) using @ yields the BdG eigen- 
value problem. The eigenfrequencies ujj come in real 
(marginally stable) or imaginary (exponentially unsta- 
ble) pairs or as complex (oscillatorily unstable) quartets. 

Our analytical approach examines perturbations of the 
time-independent GP equation ([I]) with constant poten- 
tials V(z) = Vq = and g{z) = g = ±1. The pertur- 
bations in the linear and nonlinear steps are thus Ag(z) 
and AV(z) = 'yAg(z). We introduce e = |Ag| as a small 
parameter and (to facilitate presentation) use the term 
"negative width" to describe a step with Ag < 0. When 
go = ±1, equation (JTJ has two families of (standing) soli- 
ton solutions, which are characterized by center position 
£ and chemical potential fi. The case go = — 1 yields 
bright solitons <fihs(z — £) = ?7bs sech (7/b s (z — £)), where 
r/bs = yj— 2/i and [i < 0. The case go = 1 yields dark soli- 
tons <f>ds(z - = ?7ds tanh {r) ds (z - £)), where rj ds = y/p, 
and [i > 0. 

We use a Melnikov analysis [see Supplemental Ma- 
terial (SM) for details] to determine the persistence of 
bright [3T] and dark solitary waves [35]. Stable (re- 
spectively, unstable) solitary waves exist at minima (re- 
spectively, maxima) of the effective potential M^ s . Wc 
find that bright solitary waves can, in principle, be sta- 
ble within the step. For bright solitary waves, the 
eigenfrequencies u> satisfy uj 2 = 2 i-^- M" b ' s (£o) + 0(e 2 ), 
where we compute the Melnikov function's second deriva- 
tive, M£' s (£), at the solitary- wave center £ = £o- The 
eigenfrequencies for dark solitary waves satisfy ui 2 = 
i M ds(£o) (l - if ) +0(e 2 ), where we choose the root that 
satisfies Re(«u;) > 0. In contrast to the bright solitary 
waves, stationary dark solitary waves are generically un- 
stable within the step. 

Computational Results: Existence, Stability, and Dy- 
namics. We use a fixed-point iteration scheme to iden- 
tify stationary solitary-wave solutions, and we then solve 



the BdG equations numerically to determine their cor- 
responding eigenfrequencies. We employ parameter con- 
tinuation to follow the solution branches as we vary the 
step width. 
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FIG. 2: [Color online] (Top) Maximum imaginary eigen- 
frequencies versus step width (where a negative step width 
means that Ag < 0) for (left) dark solitary waves and (right) 
bright solitary waves. We show results for the perturbation 
strengths e = 0.1 and e = 0.2. Dashed curves give results for 
analytical calculations, and solid curves give numerical calcu- 
lations using the BdG equations. (Bottom) Examples of the 
corresponding eigenfrequency spectra for e = 0.1. The inset 
in the left panel shows finite-size effects (see the main text). 

We start with the £0 = branch, which exists for all 
step widths. In Fig. [2] we show the development of the 
eigenfrequencies of this branch of solutions as a function 
of step width for both dark (left) and bright (right) soli- 
tary waves. We obtain excellent quantitative agreement 
between our results from effective-potential theory and 
those from BdG computations for the nonzero eigenfre- 
quency associated with the intrinsic (translational) dy- 
namics of the solitary wave. 

For the case of repulsive BECs (g > 0), the branch of 
solutions at £ = has a real instability for Ag < (i.e., 
AV" > 0) and an oscillatory instability for Ag > 0, and 
we capture both very accurately using effective-potential 
theory. An interesting but unphysical feature of the dark 
solitary waves is the presence of small "jumps" in the rel- 
evant eigenfrequency. These jumps are finite-size effects 
that arise from the discrete numerical approximation to 
the model's continuous spectrum [33] . 

The case of attractive BECs (g < 0) is especially in- 
teresting. A pitchfork (symmetry-breaking) bifurcation 
occurs as the step widens; it is supercritical for Ag < 
and subcritical for Ag > 0. In this case, oscillatory insta- 
bilities are not possible when translational invariance is 
broken [31] . To further probe the bifurcation, we analyze 
the phase portrait of the Newtonian dynamics [M] of the 
bright solitary wave: m off ^f = -VU(£) = 2M bs (£)/iV, 
where m e ff = 1/2. See the SM for details. 

The changes to the possible trajectories in phase space 
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suggest a viable way to investigate the bifurcation ex- 
perimentally (and hence to distinguish between narrow 
and wide steps). The presence of a step alters the 
path of a moving solitary wave, as one can see espe- 
cially by examining the wave speed. As we illustrate in 
Fig. [3j the solitary-wave dynamics depend on the num- 
ber and type of phase-plane equilibria (and hence on 
the step width). The main panel shows how one can 
use variations in v cm of a transmitted bright solitary 
wave to identify which equilibria are present. We ex- 
amine trajectories starting from the same initial condi- 
tions, (zcm(O), Ucm(O)) = (4,-0.22), for step widths of 
— 1, —1.4, and —1.8. The simplest trajectory occurs for 
the narrowest width (2zq = —1): as the solitary wave tra- 
verses the step, its speed first drops before rising again in 
the center of the step and then dropping again as the soli- 
tary wave leaves the step (due to its encounter with the 
two saddles and the center in the phase plane; see Fig. [4] 
in the SM). For wider steps, the dynamics illustrate the 
effects of the bifurcation: instead of a single peak in the 
speed, there are now two peaks separated by a well. As 
the step widens further, the two peaks move outward and 
follow the centers to the edge of the step. The maximum 
and minimum in each pair move close together in both 
v cm and t as one approaches the edge of the step. The 
solitary wave can either be transmitted (as illustrated in 
Fig. |3| or reflected by the step. (For examples of reflec- 
tion, see the SM for phase-plane trajectories not entering 
the step.) 
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FIG. 3: [Color online] (Left) Effect of the step on the move- 
ment of a bright solitary wave for three different widths for 
the GP equation and for numerical solutions of the Newtonian 
dynamics of the effective-potential (EP) equations. (Right) 
Contour plots of |^/>(.z,£)| 2 obtained by solving the GP equa- 
tion numerically for step widths of —1 and —1.8. 

Conclusions. We introduced an experimentally real- 
izable setup to study homogeneous BECs with compen- 
sating linear and nonlinear potentials. We showed that 
sound speed differences and solitary-wave formation and 
motion are useful for detecting the potential step and its 
width. We also conducted a thorough examination of the 
solitary wave stability and dynamics in this collisionally 
inhomogeneous setting. 

We found that effective-potential theory gives an excel- 
lent quantitative description of the existence and eigen- 



frequencies of both bright and dark solitary waves, and 
we used it to quantitatively track the evolution of the 
translational eigenfrequencies as a function of the step 
width. We identified a symmetry-breaking bifurcation in 
the case of attractive BECs and illustrated how the pres- 
ence of the bifurcation is revealed by the motion of soli- 
tary waves through the step region. We also found that 
stationary dark solitons are generically unstable through 
either exponential or oscillatory instabilities. 

This is a promising setup for future investigations, 
as it allows the experimentally realizable possibility of 
solitary-wave control via accurate, independent tailoring 
of external linear and nonlinear potentials. 
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Supplemental Material 

Effective- Potential Theory. We determine the persis- 
tence of a bright solitary wave by calculating when its 
center position induces its associated Melnikov function 
(i.e., perturbed energy gradient) [3T] to vanish. This 
yields the equation 



+ 



d[AV(z) 
dz 

1 d[Ag(z)] xA 



dz 



dz = 0. 



(5) 



The GP equation without a potential is spatially ho- 
mogeneous and it possesses translational and [/(l)-gauge 
symmetries; these are associated with a quartet eigenfre- 
quencies at the origin. When the former is broken (e.g., 
by the steps in V(z) and g(z)), a pair of eigenfrequencies 
leaves the origin. Their evolution makes it possible to 
examine the stability of solitary waves of the perturbed 
system. We follow these eigenfrequencies by computing 
the function 



<,(£„) = 



d 2 [AV(z)} 



dz 2 



0bs(^-?o) 



1 d 2 [Ag(z)} 

2 dz 2 



<t>t s (z-£o) 



dz . 



(6) 



which determines the concavity of the perturbed energy 
landscape and is directly associated to the eigenfrequen- 
cies of the linearization through [3T] 



uj 2 = 



(7) 



where we note that M b ' s (£ ) = 0(e). Stable (respectively, 
unstable) solitary waves exist at minima (respectively, 
maxima) of the effective potential M^ s . Hence, bright 
solitary waves can, in principle, be stable within the step. 

We compute analogous expressions for dark solitary 
waves, but the Melnikov function now needs to be renor- 
malized due to the presence of a nonzero background den- 
sity [32] • We obtain 



and 



d[AV(z)] 



dz 



fad. 



1 d[Ag(z)} 

2 dz 



[Vds ~ 4>d S { Z ~ Co)] 



dz = (8) 



d 2 [AV(z)} 



+ 



dz 2 
ld 2 [Ag(z)} 



[vis - <j>ia(z - €o)] 



dz^O. (9) 



2 dz 2 

The expression for the associated eigenfrequencies is 
1 



<S = jM^o) 1 



y 



+ 0(e 2 ), 



(10) 



where we choose the root that satisfies Re(iw) > and 
we note that M^(Co) = 0(e). 

The major difference between the spectra for dark ver- 
sus those for bright solitary waves is that the continuous 
spectrum associated with the former (due to the back- 
ground state) lacks a gap about the origin. Consequently, 
exiting along the imaginary axis is not the only way for 
eigenfrequencies to become unstable. Even when they 
exit toward the real axis, they immediately leave it as 
a result of their collision with the continuous spectrum; 
this leads to an eigenfrequency quartet. Thus, station- 
ary dark solitary waves are generically unstable within 
the step. 



Bifurcation Analysis for Bright Solitary Waves. The 
case of attractive BECs (g < 0) is especially interesting. 
A pitchfork (i.e., symmetry-breaking) bifurcation occurs 
as the step widens; it is supercritical for Ag < (i.e., 
AV" > 0) and subcritical for Ag > 0. Oscillatory insta- 
bilities are not possible when translational invariance is 
broken [5T] , 

To further probe the bifurcation, we examine the New- 
tonian dynamics |34j of the bright solitary wave: 



dt 2 



-VU(S) = 2MU0/N, 



(11) 



where m c g — 1/2. We examine phase portraits of this 
system by plotting the center-of-mass position z cm sa £ 
versus center-of-mass velocity v cm « t|. As we illustrate 
in Fig. [4j this is convenient for examining the change in 
dynamics as we alter the step width. For narrow steps 
(e.g., a width of 2zq = — 1), there is a center at z cm = 
straddling two saddle points (stars) just outside of the 
step (indicated by the dash-dotted line). When Ag < 
(i.e., AV > 0), a supercritical pitchfork bifurcation oc- 
curs at 2zq R3 —1.2, as the center at the origin transitions 
to a pair of centers separated by a saddle at the origin 
(top right panel). As the step widens further (bottom 
left panel), the heteroclinic orbit that previously enclosed 
the central three equilibria is no longer present, and the 
centers are now surrounded by homoclinic orbits ema- 
nating from the outer saddle points. Eventually, each 
outer saddle and its associated center annihilate one an- 
other (bottom right panel). When Ag > 0, the types of 
equilibria are interchanged (saddles become centers and 
vice versa). The main difference that occurs in this case 
is that solitary waves can no longer be reflected by the 
step; they are all transmitted. As one increases the mag- 
nitude of the step width from 0, there is a saddle flanked 
by two centers. At the bifurcation point, the central sad- 
dle splits into two saddles with a center between them. 
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FIG. 4: [Color online] Phase planes for Newtonian dynamics describing bright solitary waves in an attractive BEC for four 
different step widths. The thick dash-dotted lines represent the step edge. We highlight the equilibria with dots, triangles, and 
stars. The orange curves correspond to trajectories originating at equilibria. We plot example trajectories in black. 



